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ABSTRACT 

We are engaged in a programme of imaging with the STIS and NICMOS (NIC2) 
instruments aboard the Hubble Space Telescope (HST), to search for the galaxy coun- 
terparts of 18 high-redshift z > 1.75 damped Lya absorption lines and 5 Lyman-limit 
systems seen in the spectra of 16 target quasars. This paper presents the results of the 
imaging campaign with the NIC2 camera. We describe the steps followed in reducing 
the data and combining in mosaics, and the methods used for subtracting the image 
of the quasar in each field, and for constructing error frames that include the system- 
atic errors associated with the psf subtraction. To identify candidate counterparts, 
that are either compact or diffuse, we convolved the image and variance frames with 
circular top-hat filters of diameter 0.45 and 0.90 arcsec respectively, to create frames 
of summed S/N within the aperture. For each target quasar we provide catalogues 
listing positions and aperture magnitudes of all sources within a square of side 7.5" 
centred on the quasar, detected at S/N > 6. We find a total of 41 candidates of 
which three have already been confirmed spectroscopically as the counterparts. We 
provide the aperture magnitude detection limits as a function of impact parameter, 
for both detection filters, for each field. The average detection limit for compact (dif- 
fuse) sources is Hab = 25.0 (24.4) at an angular separation of 0.56" (0.79") from the 
quasar, improving to Hab = 25.5 (24.8) at large angular separations. For the brighter 
sources we have measured the half-light radius and the n parameter of the best-fit 
deconvolved Sersic-law surface-brightness profile, and the ellipticity and orientation. 

Key words: galaxies: formation - quasars: absorption lines - quasars: individ- 
ual: CS73, PC 0056-1-0125, PHL1222, PKS 0201-^113, 0216+0803, PKS 0458-02, 
PKS 0528-250, H 084H-1256, PC 0953+4749, B2 1215+33, Q 1223+1753, H 1500#13, 
Q2116-358, Q2206-1958, BR2212-1626, 2233.9+1318 



1 INTRODUCTION 

The bulk of the neutral gas in the Universe is contained in 
clouds of highest column density, Nhi > 2 x 10'^'^ cm~^, the 
damped Lya (DLA) absorbers. The incidence of these ob- 
jects with different Nhi in the spectra of background quasars 
allows one to determine the comoving density of neutral gas 
as a function of redshift or time. At z ~ 2-3 this density was 
roughly equal to the present comoving density of luminous 
stars. This observation has led to the standard interpreta- 
tion of the DLA absorbers as the gas reservoirs from which 
stars form, i.e., as the interstellar components of galaxies and 
protogalaxies (Wolfe et al. 1986) . This interpretation is sup- 
ported by global studies of chemical evolution, which relate 



the evolution of the comoving densities of stars, gas, heavy 
elements, and dust in galaxies (Pei, Fall, & Hauser 1999). 
This global approach, however, tells us nothing about the 
properties of individual DLA galaxies, such as their lumi- 
nosities, sizes, and morphologies. 

Deep imaging is needed to detect and characterize the 
stellar components of individual DLA galaxies. This would 
also provide impact parameters for absorption, which, from 
the measured incidence dn/dz, would allow the calculation 
of the space density (M0ller & Warren 1998). The combi- 
nation of information on both the stellar (emission) and in- 
terstellar (absorption) components of the DLA galaxies is 
crucial for a complete understanding of these objects. At 
present, there is much debate about the nature of the DLA 
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galaxies: on the basis of the velocity structure in metal ab- 
sorption lines, Prochaska & Wolfe (1997) argue that DLA 
absorbers are the fully-formed gaseous disks of present-day 
spiral galaxies, while Haehnelt et al. (1998) demonstrate 
that the same data can be interpreted equally well as the 
signatures of merging protogalactic clumps. Helping to re- 
solve such issues is the motivation for the imaging campaign 
described here. 

In this paper we are concerned with DLA absorbers 
at high-redshift z > 1.75. (The best imaging study of low- 
redshift DLA absorbers is by Le Brun et al 1997.) Attempts 
to image high-redshift DLA absorbers from the ground have 
mostly been unsuccessful. M0ller and Warren (1998) and 
Fynbo, M0ller, and Warren (1999) provide a summary of the 
observations and the conclusions that can be drawn from 
them. The measured impact parameters of the handful of 
detections are in the range 0.9" to 2.9", and magnitudes are 
V ^ 25. (At z = 2.5 an impact parameter of l"corresponds 
to 3.9h~^ kpc physical separation for Q,m ~ 1.0, Q,a ~ 0.0, 
and to 5.7h~^ kpc for Q,rn = 0.3, f^A = 0.7, where h = 
-ffo/100.) One may reasonably suppose that the absorbers 
that have escaped detection have similar impact parameters, 
but are fainter, or have smaller impact parameters. In either 
case the clear choice of telescope for a new imaging campaign 
is the Hubble Space Telescope, for the unrivaled depth and 
spatial resolution achievablqj. 

We are engaged in a programme of imaging with the 
STIS and NICMOS (NIC2) instruments aboard HST, to 
search for the galaxy counterparts of 18 high-redshift z > 
1.75 DLA absorption lines seen in the spectra of 16 target 
quasars. This paper presents the results of the imaging cam- 
paign with the NIC2 camera. The STIS images reach about 
2.5 AB mag. deeper than the NICMOS images, and have 
twice the resolution, and therefore are expected to be the 
most useful in terms of proportion of detections. The NIC- 
MOS observations complement the STIS observations, pro- 
viding colours, as well as luminosity profiles in the restframe 
optical. Also the NICMOS frames will be more sensitive for 
any galaxy counterparts that are particularly red. In Section 
2 we tabulate the target quasars, and the details of the ob- 
servations, and describe the steps followed in reducing the 
data and combining in mosaics. We describe the methods 
used for subtracting the image of the quasar in each field, 
and for constructing error frames that include the system- 
atic errors associated with the point spread function (psf) 
subtraction. Section 3 includes a description of the detection 
algorithm and catalogues the candidate galaxy counterparts 
found in each field, together with details of their shapes, 
sizes, and luminosity profiles, where measurable. The de- 
tection limits as a function of angular separation from the 
quasar are also provided for each field. 

In a recent preprint Kulkarni et al. (2000) present sim- 
ilar NICMOS observations of a single quasar field, where 
they detect a candidate DLA absorber counterpart at small 
impact parameter, 0.25". Their method of psf subtraction 
is different from ours. Their preprint appeared after we had 
completed this work, and the techniques were developed in- 



* Bunker et al (1999) discuss the effectiveness of blind spec- 
troscopy for detecting emission from DLA absorbers, but conclude 
that this is unlikely to be efficient even with an 8m telescope. 



dependently. We compare their method of psf subtraction 
with ours, and discuss their detection in the light of the 
detections reported here. 



2 TARGETS, OBSERVATIONS, DATA 

REDUCTION, MOSAIC COMBINATION, 
AND PSF SUBTRACTION 

2.1 Targets and observations 

The target absorbers were selected with a view to cover- 
ing a wide range of column densities in order to be able to 
investigate any correlation between column density and, for 
example, impact parameter. We wished also to cover a broad 
redshift interval, to quantify any evolution in the properties 
of the absorbing galaxies. A lower redshift limit Zabs = 1.75 
was imposed as this corresponds to the wavelength at which 
Lya enters the range observable from the ground. In fact 
many of the nearly 100 DLA absorbers so far catalogued 
lie within quite a narrow range of redshifts 1.8 < z < 2.5 
(50% of the 80 DLA absorbers listed by Wolfe et al 1995). A 
secondary consideration in selecting targets was the desire 
to avoid very bright quasars, since the faint galaxy coun- 
terparts would then be harder to detect in the glare of the 
quasar light. 

We have imaged the fields around 16 quasars. In Table 1 
we list in successive columns: 1. quasar immber — this imm- 
ber is used throughout this paper in referring to the fields 
by number, 2. quasar name — we have used the quasar name 
listed in the catalogue of Veron-Cetty and Veron (1998), 3. 
statement of whether or not the quasar has been detected at 
radio wavelengths, and 4. quasar redshift z^m- In columns 
5 and 6 are listed the column densities Nhi and redshifts 
Zabs of all the catalogued absorption lines of column density 
Nhi > lO'^^ cm~^. Column 7 lists the references, in order, 
from which the details listed in respectively columns 4, 5, 
and 6 were taken. In the final column of Table 1 each ab- 
sorber is classified as DLA if the column density is > 2 x lO^'' 
cm~^, and otherwise as LLS, standing for Lyman-limit sys- 
tem. There are 18 absorbers classified DLA, 5 classified LLS, 
and one candidate high-column density system that remains 
to be confirmed. These are the absorbers for which we aim 
to detect the counterpart galaxies. In Fig. 1 we plot column 
density against redshift for the 23 absorbers in Table 1 that 
have been confirmed with high-resolution spectroscopy. 

We used the NIC2 camera with the F160W (i.e. H band) 
filter to observe the fields around these 16 quasars, at various 
dates commencing 23 April 1998 and ending 12 Sep 1998. 
The choice of the NIC2 camera was arrived at through con- 
sideration of three factors: depth achieved, sampling of the 
psf, and field of view. For NICl the smaller pixel size car- 
ries with it a large read-out noise penalty, while the NIC3 
camera pixel size is too large to allow accurate subtraction 
of the psf. With NIC2 the psf is almost critically sampled, 
and the field of view is large enough to allow the image to 
be stepped around the array, for improved sky subtraction. 

Each field was observed for three orbits, and each orbit 
was split into two exposures of length 1280 sec, resulting 
in six exposures per quasar, totalling 7680 sec. The read- 
out mode employed was the STEP-256 MULTIACCUM se- 
quence. We used a rectangular dither pattern, with a step 
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Table 1. Redshifts and column densities of target absorption lines 
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" A strong absorption line is seen near 3480A in a low-resolution spectrum, which remains to be confirmed as Lya with a high- 
resolution spectrum (Hazard, private communication). 

The quoted absorption redshifts and column densities arc preliminary. 
Refs; [1] Bunker ct al. 1999, [2] Bunker et al., in preparation, [3] Hazard and Sargent, in preparation, [4] Lanzctta ct al 1991, [5] 
Lu and Wolfe 1994, [6] Lu ct al 1996, [7] Lu ct al 1993, [8] M0llcr ct al 1998, [9] M0llcr and Warren 1993, [10] Morton ct al 1980, 
[11] Pettini et al 1994, [12] Pettini et al 1997, [13] Veron-Cetty and Veron 1998, [14] M0ller, in preparation, [15] M0ller et al. 1994, 
[16] White et al. 1993, [17] Le Doux et al 1998. 
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size of 3", placing the quasar in the NIC2-FIX aperture for 
the first exposure of each sequence. The HST Programme 
Number of these observations is 7824. In Table 2 we hst in 
successive columns: 1. the quasar number from Table 1 (this 
is the Progammc Visit Number in the HST data archive), 2. 
quEisar name, 3. common alternatives for the quasar name 
encountered in the literature, 4. J2000 coordinates measured 
from the Digital Sky Survey (DSS) plates, 5. DSS plate used 
for the coordinates measurement, 6. The F160W Hab mag- 
nitude of the quasar measured from our data — unless stated 
otherwise all magnitudes quoted in this paper are on the AB 
system, 7. UT date at the start of the observations, 8. ref- 
erence number for locating the data in the HST archive. 
The target PC 0953-1-4749 is not visible on the DSS plate. 
To calculate the position given in Table 2 we measured the 
coordinates of a nearby galaxy on the DSS, and the offset 
to the quasar from the finding chart in Schneider, Schmidt, 
and Gunn (1991). As it happened HST was unable to ac- 
quire the guide star for this target, so the frames are badly 
trailed and could not be used. All the other observations 
were successful. 



Figure 1. Plot showing coverage of the column-density — redshift 
plane for the 23 confirmed absorbers listed in Table 1, 18 DLA 
absorption lines and 5 Lyman-limit systems. 



2.2 Data reduction, and mosaic combination 
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Table 2. Details of target quasars and observations 

Quasar Quasar Alternative Right Ascension Declination DSS Hab Date Archive 

no. name names (J2000) plate no. (total) observed rcf. 
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° Guide star acquisition failed for this target, and the images are too badly trailed to be useful. 



2.2.1 Data reduction 

We used the STSDAS task calnica, version 3.2, in the nic- 
mos package to reduce the data, in conjunction with the new 
tasks biaseq and pedsky in the nicproto package, and us- 
ing the most recent calibration frames. NICMOS data suffer 
from pedestal variations which produce offsets to the count 
level in each array quadrant. Without correction the final 
frames will bear the imprint of the flat-field frame, either 
negative or positive depending on the pedestal level. The 
tasks in the nicproto package are designed to remove this 
effect. The calnica task includes routines for subtracting 
the dark current, linearising the data, examining the mul- 
tiple reads to eliminate cosmic rays, and flat fielding the 
data. Inspection of the data in the multiple reads for an 
individual pixel is possible using the pstat task, and this 
was used to help in optimising the parameters of the differ- 
ent routines. In some cases vertical banding was visible in 
the final frames. These were removed by subtracting from 
each column the median pixel value for that column. The 
result of the first stage of the data reduction was 96 frames, 
comprising 6 exposures for each of the 16 quasars. 

The next step was to use the data frames themselves 
to remove residuals, whether multiplicative ("flat field") or 
additive ("sky subtraction"). To a first approximation the 
sky level is the same in all frames so that the residuals, 
from whatever cause, will be similar in each frame and can 
be removed by subtracting the median of the 96 frames. 
The step size between pointings, 3", was chosen to be suf- 
ficiently large that the quasar images from different frames 
would not overlap, so that objects would be successfully re- 
moved in taking the median. Subtracting the median of the 
96 frames resulted in a very noticeable reduction in the noise 
in each frame. We understand that the explanation for this 
is that the frame used for dark subtraction has substantial 
read noise, so that noise is being added to each frame at 
the stage of dark subtraction. We remove this noise by sub- 
tracting the median frame. There is a new task which pro- 
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Figure 2. Plot illustrating the variable quality of the data. The 
histogram shows the distribution of measured values of the noise 
in the sky in the 96 reduced frame, in units of DN sec~^. Taking 
into account read noise, dark noise, and Poisson noise from the 
counts in the sky, the expected noise is in the range 0.005 — 0.006 
DN sec"'^, corresponding to the peak in the histogram. As is 
evident, less than half the data lie within this range. 



duces synthetic darks of low noise, which therefore avoids 
this problem. However we were advised by the NICMOS 
group at STScI that the procedure we followed produces 
equally satisfactory results. Q 



t We investigated a second refinement to the removal of residuals, 
by considering only the six frames for an individual quasar, and 
for each frame subtracting the median of the other five frames. 
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At this stage we created a list of bad pixels by identify- 
ing discrepant high and low pixels in the median of the 96 
frames, as well as pixels where the recorded values showed 
a large dispersion in the stack of 96 frames. The bad pix- 
els were recorded as a "mask" frame for use at the stage of 
combining the six data frames from each field. 

The quality of the NICMOS data is highly variable. 
First there are a few frames which contain an excessive num- 
ber of cosmic ray hits. A related but more serious problem 
is a consequence of the passage of HST through the South 
Atlantic Anomaly (SAA) and the fact that these HgCdTe ar- 
rays suffer from remanence i.e. very high counts in a pixel are 
not completely flushed when the array is reset. Although the 
NICMOS cameras are switched off during passage through 
the SAA, the arrays experience a storm of cosmic-ray hits. In 
most cases the result of remanence of the cosmic-ray hits is 
a background with a very high noise level that affects subse- 
quent images, and declines with time. The remanence images 
of individual cosmic-ray hits are usually difficult to discern, 
presumably indicating that the whole array has been hit. To 
quantify the effect on the data we measured the standard de- 
viation in the counts in the background in each frame. Fig. 2 
plots the histogram of these noise measurements. The peak 
in the histogram marks the noise in a clean frame, unaffected 
by the SAA. As shown in Fig. 2 in 20% of our frames the 
noise is more than double that in an unadulterated frame 
and only about half the data are unaffected. 

The NICMOS pipeline produces error frames which as- 
sume a Poisson noise model. For the majority of frames this 
will be inappropriate because of the additional source of 
noise discussed above. We therefore created our own noise 
frames by adding in quadrature the Poisson contribution 
from counts above sky to the measured noise in the back- 
ground. 

2.2.2 Mosaic combination 

The pipeline task calnicb combines frames processed by 
calnica to produce a final mosaic image. We did not use 
calnicb because we wished to make a number of changes to 
the standard pipeline. The relatively poor quality of a sig- 
nificant fraction of the data demanded a more fiexible ap- 
proach that allowed each step in the combination process to 
be scrutinised. We were particularly concerned with identi- 
fying bad data, i.e. pixels with discrepant counts, either high 
or low, from whatever cause: for example due to cosmic rays 
not identified by calnica, or from transient variation of the 
quantum efficiency or dark current of any pixels. We used 
the routine driz_cr in the ditherll package (Fruchter and 
Hook 2000) to identify bad data. The first step was to align 
the six frames using integer pixel offsets, and to form the 
median, ignoring pixels in the bad-pixel mask. Each frame 
was then compared with the combination frame to identify 
discrepant data values. The bad pixels identified were com- 
bined with the original bad-pixel mask to form new individ- 
ual masks for each frame. A defect in the form of horizontal 
bars is visible in some frames. These were identified and the 
pixels added to the respective masks. Using the new masks 

This led to an increase in the noise in the background, and there- 
fore was not implemented. 



a revised combination frame was formed, and a revised iden- 
tification of bad data was made. 

The step size between pointings, 3", is by design an in- 
teger multiple of the pixel size, 0'.'075. The data frames were 
aligned to the nearest pixel, and then averaged, weighting 
by the inverse of the measured variance in the sky, and ig- 
noring masked pixels. The associated error frame was also 
created as appropriate for the weighting and mask scheme 
used. We could have chosen sub-pixel dither steps and to 
drizzle the data (Fruchter and Hook 2000). We chose not to 
because the number of pointings is so small. In the event, in 
a few frames the step pattern executed produced residuals 
from the requested 40-pixel step of up to half a pixel. The 
sub-pixel residuals therefore broaden the psf very slightly. 
Fortunately the residuals were almost identical for each tar- 
get, so the psfs of each quasar were affected in the same 
way. 

The final frames are displayed in the left hand portion 
of Figs 5 — 8 as S/N maps i.e. the image frame divided by 
the error frame. Displayed in this way the reduced S/N near 
the edges of the mosaic frames, where the number of over- 
lapping frames is smaller, is evinced by the smaller number 
of sources visible on average in these regions. Over the re- 
gions where there are only two overlapping frames the data 
are of little value. This is because the scheme used to iden- 
tify bad pixels is ineffective if the number of frames is less 
than three. In each of the pictures the bright point source 
below and to the right of the centre of the frame is the tar- 
get quasar. In a number of cases objects are visible at small 
angular separations from the quasar. These are candidate 
counterparts of the DLA absorbers. 



2.3 Psf subtraction 

To identify candidates at small angular separation from the 
quasars, and to measure their properties, we need to sub- 
tract the images of the quasars. There are three possible 
methods for psf subtraction. We could use the Tiny Tim 
software to create artificial psfs, or we could use images of 
bright stars in the archive, or we could use the 15 quasar 
images themselves. The advantages of Tiny Tim psfs over 
real psfs are that the psfs are noiseless, and can be created 
at the required x, y location on the array (the psf varies with 
X, y location), and at the precise sub-pixel position (which 
is desirable as the psf is slightly undersampled). However 
because the accuracy of the Tiny Tim psfs is limited by the 
detail in the models, real psfs are preferred if suitably placed 
high S/N images are available. ^ 

The main advantage of using images of bright stars 

i After completing this work a new version of Tiny Tim (5.0) 
was released, including modelling of the temporal variation of 
the NICMOS psf, and other improvements. We have found that 
the results obtained using the newest version of Tiny Tim are 
as good as our results obtained using the quasar images them- 
selves, as quantified by measuring the dispersion in the residuals. 
In some cases the Tiny Tim psfs are somewhat better, and in 
some cases the quasar psfs are better. As suggested in the Tiny 
Tim manual, this is probably because there are small, essentially 
random, temporal variations in the psfs, so that, by chance, one 
or the other method works better for a particular quasar. 
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is the high S/N, with the proviso that any non-linear be- 
haviour of the detector must be accurately characterised. 
The undersampling of the psf nevertheless sets a fundamen- 
tal limit to the accuracy achievable if a single star is used. 
Systematic residuals are unavoidable unless the star was po- 
sitioned at exactly the same sub-pixel position as the target. 
Other disadvantages include possible poor colour match, and 
mismatch of the psf in terms of x, y location. Primarily to 
deal with the problem of undersampling we chose to use 
the 15 quasar images themselves to subtract the psfs. Be- 
cause we have 15 targets we can use dithering to recover 
full sampling of the psf, without seriously compromising the 
5'/A'^ of the composite psf. Additional advantages of using 
the quasars are that the colours are perfectly matched (as- 
suming the range in colour is small), and that the initial x, 
y location and step pattern is the same for all the quasars. 
In addition all the observations were completed over a few 
months, so the condition of the telescope and camera was 
similar for each target. A drawback to this approach is that 
it diminishes the possibility of detecting the quasar host 
galaxies, since we will subtract the average host-galaxy pro- 
file from each image. 

The procedure followed for each quasar was to create 
a subtraction psf by averaging the images of the other 14 
quasars. Taking the relevant 14 images, the first step was 
to scale the images to the same count level, and to regis- 
ter them to the nearest half pixel. Then, using the dither 
package, the data were interlaced into pixels of side half the 
original size|^ This sub-pixelised psf was then scaled and 
registered to the quasar image by minimising x^, comput- 
ing the flux at the location of a pixel by interpolation. In 
computing the x^i pixels within 0.3" of the quasar centroid 
were excluded from the sum. Having subtracted the images 
of all 15 quasars, mask frames were created by adding to 
the mask any galaxies found close to the quasar, as well as 
any remaining bad pixels visible. The whole procedure was 
then iterated, using the masks to exclude these pixels in cre- 
ating the psfs, and in making the fits. In this second stage 
the quasar PHL 1222 was excluding from the averaging in 
creating subtraction psfs because of excess emission clearly 
visible. This excess emission is centred on the quasar and is 
likely to be from the host galaxy of the quasar. 

Sub-pixelisation of the psf has two effects. First it im- 
proves the sampling and so reduces systematic errors in sub- 
traction near the centre of the psf. However the resulting psf 
is noisier. Because of the limited number of quasars used the 
psf subtraction adds noise to the frames. This is more no- 
ticeable for the brightest quasars, since the subtraction psf is 
formed by scaling up the images of fainter quasars. The pro- 
cedure could be improved, therefore, by using a sub-pixelised 
psf for subtraction near the centre of the quasar image, but 
using full-sized pixels away from the centre where the issue 
of sampling is less important. 

The approach Kulkarni et al. (2000) followed was to use 
single bright stars for psf subtraction. They provide an ex- 
tensive and detailed analysis of the sources of uncertainty 
using this approach. In section 3.2.3 we compare the accu- 
racy of our psf subtraction with that of Kulkarni et al. 

§ In dither parlance interlacing corresponds to a delta-function 
drizzle footprint. Any finite-sized footprint broadens the psf 



2.4 Calculation of error frames 

The final stage before searching for candidates was to com- 
pute revised la error frames for each quasar. The purpose 
of this is to be able to measure the significance of candidate 
detections. The revised error frames must include the con- 
tribution of random and systematic errors associated with 
the psf subtraction. We first computed S/N frames by di- 
viding the subtracted frames by the current error frames. 
Ideally, if the subtraction psf were both noiseless, and per- 
fectly accurate, then the distribution of counts in the S/N 
frames, outside remaining objects, would be Gaussian with 
standard deviation cts/jv = 1-0. Because the subtraction psf 
is not noiseless, to account for the added noise we scaled the 
error frame by the measured a^/jv in the sky, and then cre- 
ated new S/N frames. Even then (Js/ni although now unity 
in the sky, increased towards the centre of the (subtracted) 
quasar. The increase was greater for the brighter quasars. 
The reason for this can be understood as follows. For bright 
objects, such as the quasars in this programme, the errors 
are dominated by object photon noise i.e. \/iV, where A'' is 
the number of detected photons. Suppose that the psf has 
a fractional accuracy / i.e. the typical residual for any pixel 
after psf subtraction is fN. Then in the S/N frame the resid- 
uals will be visible as deviations of S/N — fVN, i.e. larger 
where the counts are larger, and larger for brighter quasars. 

To quantify these systematic errors, for each quasar we 
measured as/N in annuli centred on the quasar, clipping out 
real objects. To these results we fit a function of the form 
""s/iv = max(1.0 -|- e~^^^{a — 6Hab), 1.0). The error frames 
were then rescaled by the fitted function. In this equation 
r is the distance from the quasar centroid in pixels, R is 
a scale length, a and b are constants, and Hab is the ap- 
parent magnitude of the quasar. The fit of the function is 
illustrated in two plots in Fig. 3. These plots indicate that 
for quasars fainter than Hab = 18.5 it is possible to achieve 
psf subtraction with negligible systematic error, but that the 
systematic errors increase rapidly at brighter magnitudes. 

It is noticeable from the right hand plot in Fig. 3 that 
the scatter increases at smaller radii, so that the scaling 
becomes less reliable near the centre. From an inspection of 
the residuals for the different targets we judged that inside a 
radius of 4.5 pixels (0.34") the psf subtraction is unreliable, 
and we set the data to zero inside this radius. 



3 DETECTION OF CANDIDATES, AND 
CATALOGUES 

3.1 Detection 

In this subsection we describe the algorithm used to pro- 
duce a list of candidates. We restricted the search to a box 
of side 7.5" centred on the quasar. A commonly used pro- 
cedure for detecting sources is to convolve the frame with 
a filter, typically a Gaussian, and to identify peaks in the 
convolved frame. There were three reasons why we did not 
follow this approach. The first is that the noise in the sub- 
tracted frames is larger near the quasar centre, so we need 
a detection threshold that varies as the noise varies. The 
second is that objects with surface-brightness profiles very 
different from the filter profile will not be detected. So, for 
example, low-surface brightness galaxies could be missed. 
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Figure 3. Illustration of tlie rescaling of the Poisson errors to account for systematic errors associated with the psf subtraction. At each 
radius for each quasar the standard deviation ag/ff in the S/N frame is measured. A function of the form 1.0 + e~''/^(a — 6Hab) is 
fit to the data. To plot all the data at once we apply the results of the fit to rescale the measured values to a specific radius r = 5 
pixels (0.375") (left-hand plot), and a specific magnitude Hab = 18.0 (right-hand plot). The left-hand plot then shows how the scaling 
increases for brighter quasars, and the right-hand plot shows how the scaling increases at smaller radii. The solid lines axe the fitted 
function. 



The third is that wc wish to quote detection hmits in a way 
such that it is possible easily to compute whether a hypo- 
thetical galaxy (of arbitrary shape, profile, and magnitude) 
would have been detected. 

Our method is to convolve with detection filters in a 
manner that produces complete samples of objects that have 
S/N above a specified threshold, as measured within circular 
apertures of specified sizes. The data and variance frames are 
convolved with circular top-hat filters, of size equal to the 
requested aperture. The frame formed by taking the ratio 
of the convolved data frame and the square root of the con- 
volved variance frame contains, at each pixel, the summed 
S/N within the aperture of any source centred at that pixel. 
This ratio frame is then searched for peaks above a chosen 
S/N threshold. 

We first ran a large-box median filter over the frames, 
and subtracted this off, to remove any remaining gradual 
variations in the background. We used two detection aper- 
tures, of diameter 6 pixels (0.45") and 12 pixels (0.90"). The 
smaller aperture was designed to detect compact sources, 
and the larger aperture was designed to find objects of lower 
surface-brightness that had been missed by the first filter. 
The S/N detection threshold was selected from a considera- 
tion of the distribution of the depths of negative peaks in the 
aperture S/N frames. With the exception of one field there 
were no negative detections in any frame below S/N = — 6 
for either detection filter. For the quasar Q2116 — 358 the 
psf subtraction was less successful than for the other tar- 
gets, and there are negative residuals at the position of two 
of the diffraction spikes with S/N ~ —8. With this in mind 



we set our detection threshold at S/N = 6 with the expec- 
tation that nearly all the objects detected will be real, with 
the caveat that sources located on diffraction spikes will be 
less reliable. In fact with the exception of very red objects 
we expect that all real objects in our candidate list will be 
detected at higher S/N in our STIS data. Therefore we will 
reassess the reliability of the candidate list on completion of 
the list of candidates detected in the STIS frames. 

Inner radius limit. Because the psf subtraction inside 
r = 0.34" was judged unreliable (section 2.4) the inner ra- 
dius to which objects can be detected is equal to the sum of 
0.34" and the aperture radius. Therefore for the small aper- 
ture the inner radius limit is 0.56", and for the large aperture 
the inner radius limit is 0.79". 

Detection limits. The candidate list, then, was con- 
structed by identifying any peaks higher than 6.0 in the 
aperture S/N frames. Real objects in the frames will be 
detected at higher S/N with the filter that most closely 
matches the galaxy surface brightness profile. Therefore if 
a peak was found at the same location for both filters we 
eliminated the detection of lower S/N, and classified the 
object as compact (C) or diffuse (D) for detection with the 
small or large aperture filter respectively. At the location of 
each peak wc integrated the counts within the aperture, and 
transformed this to an aperture Hab magnitude using the 
photometry zero point provided in the file header (the data 
were calibrated by converting the count rates (DN/scc) to 
Jy by mulitplying by the factor 2.070 x lO""). Then for each 
frame and for both filters we measured the average detec- 
tion limit in annuli at several different radii. These values 
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Figure 4. Plot of relation between Hab aperture magnitude and impact parameter for the 41 candidate DLA galaxies, taken from Table 
4. Solid symbols identify compact objects (C in Table 4), and open symbols diffuse objects (D in Table 4): circles mark candidates, 
triangles mark sources confirmed as the DLA counterpart. The curves mark the detection limits averaged for the 15 fields, solid for 
compact objects and dashed for diffuse sources. The vertical solid and dashed lines mark the inner radius limit used for the search for 
objects, for the small aperture (i.e. compact objects) and large aperture (i.e.diffuse objects) respectively. Candidates with centroids inside 
these limits are not considered reliable and have been excluded from the catalogue in Table 4. 



are listed In Table 3. In Fig. 4 we plot, for both filters, the ra- 
dial profile of the detection limits averaged for the 15 fields. 
The average detection limit for compact (diffuse) sources is 
Hab = 25.0 (24.4) at an angular separation of 0.56" (0.79") 
from the quasar, improving to Hab = 25.5 (24.8) at large 
angular separations. 

Aperture magnitudes. All the magnitudes of the candi- 
dates (section 3.2) are quoted as aperture magnitudes i.e. 
we have simply summed the counts within the aperture and 
no aperture correction has been applied (aperture correc- 
tions are, nevertheless, discussed in section 3.2.1). We have 
not attempted to extrapolate to total magnitudes, since this 
would be unreliable for the faint sources. The use of aper- 
ture magnitudes makes it very simple to compute whether 
or not a hypothetical galaxy of specified surface brightness 
profile and impact parameter would have been detected in a 
particular field. One first computes the aperture magnitude 
of the galaxy, for the small and large apertures. For the hy- 
pothesised impact parameter one compares the two aperture 
magnitudes against the detection limits provided in Table 3 
for the field in question. If the galaxy were brighter than 
either detection limit it would have been included in our 
catalogue. 



3.2 Catalogue 

The catalogue of detections is provided in Table 4. For each 
quasar field the candidates are listed in order of angular sep- 
aration from the quasax. Listed in the first three columns of 
Table 4 are the quasar number, quasar name, and the candi- 
date number. The candidate numbering scheme is explained 
as follows. Taking as example candidate N-15-3C, the 'N' 
stands for NICMOS (we will provide a similar list 'S' for 
the STIS images), 15 is the quasar number, the 3 indicates 
that it is the third nearest candidate to the quasar in that 
field, and the C states that it is compact. We find a total 
of 41 candidates. The magnitudes and impact parameters 
are compared against the average detection limits in Fig. 4. 
The positions of the candidates arc illustrated in the finding 
charts in the right-hand panels in Figs 5 — 8. The remaining 
columns in Table 4 list successively the coordinates, aper- 
ture magnitude, detection S/N, angular separation from the 
quasar, and the position angle from the quasar. The NIC- 
MOS data header files contain an astromctric solution that 
provides coordinates that are accurate in a relative sense, for 
the small angular separations with which we are concerned. 
The coordinates were shifted to absolute values by adopting 
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Table 3. Aperture magnitude Hab average detection limits as a function of impact parameter for each field, for compact and diffuse 
sources. 

No. Quasar compact or impact parameter (arcsec) 

diffuse 0.3 0.4 0.5 0.6 0.8 1.0 1.5 2.0 3.0 

1 CS73 C 23.77 24.34 24.94 25.21 25.44 25.58 25.71 25.72 25.74 

D 22.65 22.75 23.47 24.02 24.63 24.80 24.95 24.98 25.00 

2 PC 0056+0125 C 23.89 24.44 24.95 25.16 25.37 25.48 25.59 25.62 25.63 

D 22.80 22.92 23.60 24.10 24.57 24.72 24.84 24.88 24.89 

3 PHL1222 C 23.51 24.08 24.74 25.08 25.39 25.59 25.75 25.78 25.79 

D 22.35 22.45 23.20 23.84 24.55 24.79 25.00 25.04 25.05 

4 PKS 0201+113 C 25.33 25.56 25.71 25.76 25.78 25.80 25.81 25.81 25.81 

D 24.46 24.54 24.82 24.95 25.03 25.05 25.07 25.07 25.07 

5 0216+0803 C 23.42 24.02 24.71 25.03 25.31 25.49 25.64 25.69 25.70 

D 22.23 22.34 23.13 23.75 24.48 24.71 24.89 24.94 24.96 

6 PKS 0458-02 C 25.24 25.39 25.51 25.54 25.55 25.56 25.57 25.57 25.57 

D 24.36 24.45 24.66 24.75 24.80 24.80 24.82 24.83 24.83 

7 PKS 0528-250 C 23.13 23.70 24.33 24.65 24.95 25.16 25.33 25.39 25.40 

D 21.97 22.07 22.81 23.41 24.12 24.37 24.59 24.64 24.66 

8 H 0841+1256 C 23.81 24.20 24.61 24.78 24.98 25.10 25.21 25.24 25.25 

D 22.75 22.86 23.45 23.84 24.19 24.34 24.46 24.50 24.51 

10 B2 1215+33 C 24.55 24.97 25.31 25.45 25.57 25.64 25.68 25.70 25.70 

D 23.48 23.58 24.16 24.52 24.80 24.88 24.93 24.95 24.96 

11 Q 1223+1753 C 23.25 23.77 24.29 24.55 24.82 25.00 25.16 25.20 25.21 

D 22.15 22.26 22.94 23.47 24.00 24.21 24.41 24.46 24.47 

12 H1500#13 C 25.27 25.49 25.65 25.69 25.70 25.71 25.72 25.72 25.72 

D 24.41 24.48 24.75 24.88 24.95 24.97 24.98 24.98 24.98 

13 Q 2116-358 C 22.72 23.27 23.86 24.22 24.57 24.82 25.05 25.11 25.12 

D 21.46 21.56 22.40 22.97 23.72 24.01 24.31 24.36 24.38 

14 Q 2206-1958 C 22.49 23.06 23.92 24.34 24.73 25.04 25.30 25.38 25.40 

D 21.23 21.34 22.10 22.89 23.87 24.21 24.55 24.63 24.65 

15 BR 2212-1626 C 23.65 24.11 24.58 24.76 24.96 25.09 25.22 25.25 25.25 

D 22.61 22.71 23.36 23.77 24.17 24.33 24.46 24.50 24.51 

16 2233.9+1318 C 24.02 24.56 25.14 25.39 25.60 25.73 25.83 25.85 25.86 

D 22.88 22.99 23.68 24.26 24.81 24.95 25.08 25.11 25.12 



for the coordinates of the quasar the value measured from 
the DSS plate, listed in Table 2. 



3.2.1 Surface-brightness profiles 

For the brighter candidates wc have measured the decon- 
volved surface brightness profile using the techniques de- 
scribed by Warren et al (1996). We use the Sersic model 
where the surface brightness as a function of radius r is 
E = Se exp{— B(n)[(r/re)^^" — 1]}. The parameter n char- 
acterises the shape of the profile: n = 1 is the exponential 
profile and n = 4 the de Vaucouleurs profile. This param- 
eterisation is particularly useful, therefore, as the value of 
n can be used to classify faint galaxies into "early" and 
"late" types. The parameter is the half-light radius, Ee 
is the surface brightness at the half-light radius, and B{n) 
is a constant for particular n. Ciotti and Bertin (1999) 
provide the series asymptotic solution for B{n), and we 
used the approximation provided by the first four terms 
B{n) = 2n - 1/3 + 4/(405n) + 46/(25515n^) which is ac- 
curate to better than one part in 10® over the range of n of 
interest. 

We created a psf using the Tiny Tim software, with 
pixels one-third the size of the NIC2 pixels, to ensure ad- 
equate sampling. The fitting proceeds by creating a two- 
dimensional galaxy surface-brightness profile with the same 
small pixel size, convolving with the psf, rebinning to the 



full pixel size and computing the of the fit. The best fit 
is found by minimisation on the seven parameters x, y, 
Eg, Te, n, orientation, and ellipticity. In Table 5 we provide 
the details of the fits to all the candidates in Table 4 with 
S/N > 18. Also listed there arc the total magnitudes, com- 
puted by extrapolating the models to infinite radius, and 
the aperture corrections i.e. the difference between the to- 
tal and aperture magnitudes. It is of interest to compare 
the aperture corrections for these galaxies with the values 
for a point source. For a point source the aperture correc- 
tion for the small aperture (used for compact objects) is 
0.444 mag. Excluding the two point sources listed in Table 
5, there are seven compact objects, for which the mean aper- 
ture correction is 0.66 mag., with a scatter of 0.16 mag. For 
a point source the aperture correction for the large aperture 
(used for diffuse objects) is 0.134 mag. There are four dif- 
fuse objects listed in Table 5, for which the mean aperture 
correction is 0.60 ± 0.17 mag. 



3.2.2 Results 

We now comment on the results for the individual fields in 
turn. 

1 . CS 73, Fig. 5. The spectrum of this source shows a 
LLS at = 1.8862 and a DLA at ^; = 2.0713 (Table 1), but 
no candidate counterparts were found in this field. 
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Table 4. Coordinates and aperture magnitudes of candidate DLA galaxies 
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21 
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21 
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14 


Q 2206-1958 


N-14-1C 


22 


08 


52.1019 


-19 


43 


58.456 


25.11 


6.5 


1.13 


23.6 


Notes 1, 3 






N-14-2C 


22 


08 


52.0480 


-19 


43 


58.192 


24.53 


12.1 


1.33 


-13.2 








N-14-3C 


22 


08 


52.1966 


-19 


43 


58.213 


25.27 


7.2 


2.20 


54.3 








N-14-4C 


22 


08 


51.8986 


-19 


43 


59.055 


25.11 


7.8 


2.44 


-79.9 








N-14-5D 


22 


08 


52.2106 


-19 


44 


02.484 


24.64 


6.3 


3.54 


146.3 




15 


BR2212-1626 


N-15-1D 


22 


15 


27.3782 


-16 


11 


32.611 


24.31 


6.7 


1.14 


100.4 








N-15-2C 


22 


15 


27.1966 


-16 


11 


32.784 


21.22 


207.5 


1.53 


-104.7 








N-15-3C 


22 


15 


27.3681 


-16 


11 


34.450 


22.99 


47.5 


2.25 


154.6 








N-15-4C 


22 


15 


27.5349 


-16 


11 


32.041 


25.14 


7.0 


3.38 


83.7 








N-15-5C 


22 


15 


27.1766 


-16 


11 


29.211 


25.37 


6.2 


3.61 


-29.0 




16 


2233.9+1318 


N-16-1D 


22 


36 


19.2800 


13 


26 


18.405 


25.05 


6.6 


2.78 


158.6 


Note 4 






N-16-2C 


22 


36 


19.4239 


13 


26 


22.236 


25.12 


12.2 


3.32 


68.3 





1. The three sources N-5-1C, N-13-1D, and N-14-1C lie on diffraction spikes and should therefore be considered less certain/reliable 
than the other sources. 

2. This source is the DLA at z = 2.8110 confirmed by M0ller and Warren (1993), and so has been included even though it has 

S/N < 6. 

3. This source has been confirmed by us as the DLA at z = 1.9205 (M0ller et al, in preparation). 

4. This source is the LLS at z = 3.1501 confirmed by Djorgovski et al. (1996). 
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Figure 5. N.B. This is not the full resolution figure which is available at astro.ic.ac.uk/Research/extragal/dla.html. 

Reduced frames showing the images of four of the targets before and after psf subtraction. The pixel size is 0'.'075. The targets are, 
from top to bottom, nos 1. CS 73, 2. PC 0056+0125, 3. PHL 1222, 4. PKS 0201+113. In each case the left image is the pixel S/N frame 
formed by dividing the data frame by the error frame, and shows the full mosaic field 25'.'0 X 22'.'2. The compass points N with the 
arm to the E. The length of the arm is 3". The central image is a x3 zoom of the 7'.'5 x 7'.'5 region centred on the target, after psf 
subtraction. The vertical bar is of length 1". This shows the S/N frame for the small detection filter i.e. the S/N within an aperture of 
diameter 0.45"centred at that point. Pealcs above S/N = 6 are compact candidates. To detect diffuse candidates a corresponding frame 
for an aperture of diameter 0.9" was constructed. The right-hand frame identifies the candidates, listed in Table 4. The symbol size scales 
with brightness. Filled symbols are compact candidates, open symbols diffuse. Brighter objects arc plotted at the correct deconvolved 
orientation and ellipticity (Table 5), while fainter sources are plotted as circles. The cross marks the quasar position. 
© 0000 RAS, MNRAS 000, 000-000 
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Figure 6. N.B. This is not the full resolution figure which is available at astro.ic.ac.uk/Research/extragal/dla.html. 

Reduced frames showing the images of four of the targets before and after psf subtraction. The pixel size is 0'.'075. The targets are, 
from top to bottom, nos 5. 0216+0803, 6. PKS 0458-02, 7. PKS 0528-250, 8. H 0841+1256. hi each case the left image is tiic pixel S/N 
frame formed by dividing the data frame by the error frame, and shows the full mosaic field 25'.'0 X 22'.'2. The comijass points N with 
the arm to the E. The length of the arm is 3". The central image is a x3 zoom of the 7'.'5 x 7'.'5 region centred on the target, after psf 
subtraction. The vertical bar is of length 1". This shows the S/N frame for the small detection filter i.e. the S/N within an aperture of 
diameter 0.45"centred at that point. Pealcs above S/N = 6 are compact candidates. To detect diffuse candidates a corresponding frame 
for an aperture of diameter 0.9" was constructed. The right-hand frame identifies the candidates, listed in Table 4. The symbol size scales 
with brightness. Filled symbols are compact candidates, open symbols diffuse. Brighter objects are plotted at the correct deconvolved 
orientation and ellipticity (Table 5), while fainter sources are plotted as circles. The cross marks the quasar position. 
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Figure 7. N.B. This is not the full resolution figure which is available at astro.ic.ac.uk/Research/extragal/dla.html. 

Reduced frames showing the images of four of the targets before and after psf subtraction. The pixel size is 0'.'075. The targets are, 
from top to bottom, nos 10. B2 1215+33, 11. Q 1223+1753, 12. H1500#13, 13. Q 2116-358. In each case the left image is tlic pixel S/N 
frame formed by dividing tlie data frame by the error frame, and shows the full mosaic field 25'.'0 X 22'.'2. The comijass points N witli 
the arm to the E. The length of the arm is 3". The central image is a x3 zoom of the 7'.'5 x 7'.'5 region centred on the target, after psf 
subtraction. The vertical bar is of length 1". This shows the S/N frame for the small detection filter i.e. the S/N within an aperture of 
diameter 0.45 "centred at that point. Peaks above S/N = 6 are compact candidates. To detect diffuse candidates a corresponding frame 
for an aperture of diameter 0.9" was constructed. The right-hand frame identifies the candidates, listed in Table 4. The symbol size scales 
with brightness. Filled symbols are compact candidates, open symbols diffuse. Brighter objects are plotted yat the correct deconvolved 
orientation and ellipticity (Table 5), while fainter sources are plotted as circles. The cross marks the quasar position. 
© 0000 RAS, MNRAS 000, 000-000 
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Table 5. Sersic profile fits to brighter candidates 
No. Quasar Candidate 



orient 11 

■•' (E of N) 



ellii) 



Had 

(loial) 



AH VP, coiniiiciits 
(ap. (-iirr.) 



2 




N-2-4D 


0.212" 


Q 
O 




IVT Q 1 


n 1 1 o' 






N-3-2C 








N-3-3D 


0.437; 


4 


PKS 0201+113 


N-4-1C 


0.203; 


6 


PKS 0458-02 


N-6-5C 


0.103; 


8 


H 0841+1256 


N-8-1C 


0.068; 






N-8-3D 


0.37i; 


10 


B2 1215+33 


N-10-2D 


0.256; 


13 


Q 2116-358 


N-13-2C 


0.130; 






N-13-3C 


0.170; 


15 


BR2212-1626 


N-15-2C 








N-15-3C 


0.115; 



-0.008 
-0.006 
-0.006 



-0.035 
+0.011 
-0.009 
-1-0.008 
-O.OOS 
-1-0.005 
-0.005 
-1-0.018 
-0.015 
-1-0.035 
-0.024 
+0.026 
-0.016 
+0.010 
-0.009 



-0.006 
-0.006 



2.57; 
0.65 

0.84 
1.21 
0.41 
2.69 
0.61 
0.93 



-1-0.15 
0.14 
+0.17 
0.16 

-1-0.21 
-0.16 
-1-0.15 
0.13 
-10.25 
0.19 
-10.50 
0.40 
-10.12 
-0.10 
-10.36 
-0.25 
-1-0.82 



^•^•^-0.54 
n fjr;+0.20 



1.28 



-10.22 



64 4+°-^ 
-27.il;;; 

-1.9+^-" 



3.0 

4.5 
4.5 
+ 1.7 
-1.7 
+5.1 
-4.9 
+2.5 
-2.6 
-13.4 
-12.8 

29.7+8-2 



53.0 
-53.1 
-78.8 

72.4 
-86.0 



89.3 



+2.9 
2.9 



0.70 
0.05 

0.21 
0.38 



+0.01 
0.01 
+0.09 



+0.07 
-0.07 
+0.03 
0.03 

^•°^-0.07 

70+° °3 

'-'•"-'-0.03 

'-'•'^^-0.05 
Q 73+0.04 
''-'-0.04 
S+0.09 
-0.10 
S+0.07 
-0.08 



0.28: 
0.28"* 



0.50 



+0.04 



21.06 
23.35 

22.54 
22.40 
23.79 
21.95 
22.11 
23.21 
22.74 
22.74 

22.37 



0.45 
0.51 

0.84 
0.96 
0.50 
0.57 
0.62 
0.50 
0.72 
0.77 

0.62 



Note 1 

point source (Q 0151+0448B) 
faint companions masked 



point source 



1. This source is almost round. In consequence the confidence limits on the orientation are indeterminate, as is the lower confidence 
limit for the ellipticity. 



2. PC 0056+0125, Fig. 5. There is a strong DLA line 
at z = 2.7771 in the spectrum of this quasar. Bearing in 
mind the small impact parameters of the few confirmed 
DLA counterparts we consider the nearest candidate, N-2- 
IC which lies at an angular separation of 2.32", to be the 
most likely of the four candidate counterparts listed in Ta- 
ble 4. The large ellipticity, 0.70, of the candidate N-2-4D 
suggests that this is a late-type galaxy viewed at a high in- 
clination angle. The orientation of the image is only 10° from 
the line joining the quasar to the galaxy. Therefore a possi- 
ble alternative explanation of the DLA line is the presence 
of an extended gaseous disk around this galaxy. 

3. PHL 1222, Fig. 5. The redshift of the DLA z = 
1.9342 is close to the redshift of the quasar z = 1.922. 
The DLA absorbs much of the quasar Lya emission line, 
and to measure the column density it was necessary to 
model the unabsorbed Lya emission line profile (M0llor ct 
al 1998). There are significant residuals visible after the 
psf subtraction, which are approximately circularly sym- 
metric, which wo suggest arc most likely to bo duo to the 
quasar host galaxy. The DLA absorber has been detected in 
Lya emission by Fynbo ot al. as an extended object 6 arc- 
sec across, centred 0.9"E of the quasar. These authors sug- 
gest that the source of ionising photons is the nearby quasar 
Q0151+0448B, rather than young stars. This quasar lies at 
an angular separation of 3.32", and is listed in Table 4 as 
source N-3-2C. We do not find a source coincident with the 
centre of the Lya emission, consistent with the photoionisar 
tion explanation. 

4. PKS 0201 + 113, Fig. 5. The spectrum of this source 
shows a DLA absorber at « = 3.3875. There is a rela- 
tively bright compact source N-4-1C with Hab = 23.36 only 
1.0l"from the quasar, which is a good candidate for the 
counterpart. This object has n — 1.21 i.e. the surface- 
brightness profile approximates an exponential. 

5. 0216+0803, Fig. 6. The spectrum of this source 
shows a confirmed DLA absorber at z = 2.2930. There is 
also a high-column density system s± z = 1.7688 confirmed 
as such by Lu et al. (1996) on the basis of the detection 
of the weak metal lines of Fell A2260, Mnll A2576, and 
Nill A1751. The value for the HI column density however 



comes from the the Lya line equivalent width measured by 

Lanzetta et al (1991) from a low resolution spectrum, and is 
therefore uncertain. The nearest candidate counterpart N-5- 
IC lies at an angular separation of 1.41". As noted in Table 4 
this object lies on a diffraction spike, and therefore should be 
considered less reliable than other sources, especially given 
the relatively low S/N of 6.3. 

6. PKS 0458-02, Fig. 6. This quasar shows extended 
radio emission. Analysis of the 21 cm line in absorption al- 
lowed Briggs et al (1989) to place a lower limit on the size 
of the DLA absorber {z = 2.0395) of 8/i"^ kpc. The near- 
est candidate counterpart is the faint Hab = 24.60 diffuse 
object N-6-1D, at an angular separation of only 0.86 arcsec. 
This candidate is the closest to the quasar of all the candi- 
dates listed in Table 4, and very close to our inner radius 
limit of 0.79 arcsec for the detection of diffuse objects. The 
object docs not stand out clearly in Fig. 6 (middle column) 
as this is the detection frame for compact objects. We see no 
reason to doubt the reality of this object. The source N-6-3C 
lies on the diffraction spike of a nearby bright star, and the 
measured magnitude will overestimate the brightness of this 
source. 

7. PKS 0528-250, Fig. 6. There are two DLA absorbers 
in the spectrum of this quasar at z = 2.1408 and z = 2.8110. 
The galaxy counterpart of the higher redshift system was 
discovered by M0ller and Warren (1993) and has been ex- 
tensively studied (Warren and M0ller 1996, M0ller and War- 
ren 1998). This source is detected in our NICMOS image at 
S/N = 5.0, and we have included the object in Table 4 
even though it is fainter than the detection limit S/N = 6. 
The STIS image of the source is presented in M0ller et al 
(in preparation). The z = 2.8110 DLA system absorbs the 
quasar Lya emission line. The column density quoted for this 
absorber in Table 1 is taken from MoUcr and Warren (1993) 
who attempted to account for the quasar Lya emission in 
making the fit (in the same manner as with PHL 1222). 

8. H 0841+1256, Fig. 6. The redshift of H 0841+1256 is 
uncertain. The object was discovered by Hazard in a search 
of an objective prism plate, and may be a BL Lac as no emis- 
sion lines have been detected. The redshift is estimated from 
the wavelength of the onset of absorption in the Lya forest. 
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Figure 8. N.B. This is not the full resolution figure which is available at astro.ic.ac.uk/Research/extragal/dla.html. 

Reduced frames showing the images of three of the targets before and after psf subtraction. The pixel size is 0'.'075. The targets are, from 
top to bottom, nos 14. Q 2206-1958, 15. BR2212-1626, 16. 2233.9+1318. In each case the left image is the pixel S/N frame formed by 
dividing the data frame by the error frame, and shows the full mosaic field 25'.'0 x 22'.'2. The compass points N with the arm to the E. The 
length of the arm is 3". The central image is a x3 zoom of the 7'.'5 x 7'.'5 region centred on the target, after psf subtraction. The vertical 
bar is of length 1". This shows the S/N frame for the small detection filter i.e. the S/N within an aperture of diameter 0.45"centred at 
that point. Peaks above S/N = 6 are compact candidates. To detect diffuse candidates a corresponding frame for an aperture of diameter 
0.9"was constructed. The right-hand frame identifies ttie candidates, listed in Table 4. The symbol size scales with brightness. Filled 
symbols are compact candidates, open symbols diffuse. Brighter objects are plotted at the correct deconvolved orientation and ellipticity 
(Table 5), while fainter sources are plotted as circles. The cross marks the quasar position. 



but this cannot be measured accurately as it coincides with 
the wavelength of the z = 2.4764 DLA absorber, which also 
absorbs any Lya emission that might be present. There are 
two confirmed DLA absorbers in the spectrum of this object, 
at 2 = 2.3745 and z = 2.4764. A candidate DLA absorption 
line near 3480A remains to be confirmed. 

The source N-8-1C was first detected by Aragon- 
Salamanca, Ellis, and O'Brien (1996) in a K-band image ob- 
tained at the United Kingdom Infrared Telescope (UKIRT). 



They measured K = 19.9 (Johnson) for this source. Taking 

the total magnitude from Table 5 and converting to the 
Johnson system, our measurement gives H = 20.7, and a 
colour H - K = 0.8. 



9. PC 0953+4749. As stated earlier the NICMOS ob- 
servations of this source failed. The redshifts and column 
densities of the three DLA absorbers were measured from 
a spectrum taken with the Keck I telescope. The spectrum 
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will be presented in Bunker et al. (in preparation) together 
with WFPC2 observations of the field. 

10. B2 1215+33, Fig. 7. There is a single DLA absorber 
at z = 1.9990 in the spectrum of this quasar. This field 
was also observed by Aragon-Salamanca ot al (1996) who 
report the detection of a candidate with K = 20.1 (John- 
son) located 1.3" E of the quasar. This candidate therefore 
is 0.2 mag fainter in K than N-8-1C. However N-8-1C has 
an aperture magnitude of Hab = 22.5 whereas we do not de- 
tect their fainter candidate above our detection thresholds of 
Hab = 25.7 and 24.9 for compact and diffuse sources respec- 
tively. If the object had a similar colour and surface bright- 
ness profile to N-8-1C wc would have expected to find a 
compact source with an aperture magnitude of Hab = 22.7. 
This implies that their candidate, if real, is exceptionally 
red. 

11. Q 1223+1753, Fig. 7. There is a single DLA ab- 
sorber at z = 2.4658 in the spectrum of this quasar, but 
no candidate counterparts have been detected. Several of 
the images of this field suffer from high background noise 
(§2.2.1). 

12. H1500#13, Fig. 7. There is a DLA absorber at 
z — 3.1714 seen in the spectrum of this source. After psf 
subtraction a strong residual to the SW and close to the 
quasar remains, which is clearly visible in Fig. 7. The cen- 
troid of the residuals lies inside our detection radius limit, 
and the tangential orientation of the residual leads us to 
believe that it is not duo to a galaxy counterpart of the ab- 
sorber, but could possibly be emission from the quasar host 
galaxy. 

13. Q2116-358, Fig. 7. The spectrum of this quasar 
shows a LLS absorption line at z = 1.9966. This is the sec- 
ond brightest quasar in our list. The psf subtraction for this 
quasar is particularly poor, leading to strong negative resid- 
uals for two of the diffraction spikes. This is the only object 
for which such a pattern of residuals is seen, suggesting that 
the observing configuration was difi'erent in some way. Al- 
though we have been unable to identify a cause we note that 
this quasar was the first to be observed in our campaign. A 
ring of excess emission is seen in Fig. 7 which could be due to 
the host galaxy of the quasar. The candidate source closest 
to the quasar lies on a difi^raction spike and therefore should 
be considered less reliable than other sources. 

14. Q2206-1958, Fig. 8. The nearest source to the 
quasar N-14-1C, which lies at an offset of 1.13", has been 
confirmed by us as the counterpart to the DLA absorber at 
z = 1.9205 (M0ller et al., in preparation). The source lies 
on a diffraction spike. The psf subtraction for this source is 
comparatively poor, so the uncertainty on the magnitude is 
larger than that indicated by the S/N quoted in Table 4. 
The spectrum of this quasar shows a second DLA absorber 
at 2 = 2.0762. 

15. BR 2212- 1626, Fig. 8. There is a LLS absorber at 
z — 3.6617 in the spectrum of this quasar. There are two 
relatively bright candidates in this field. With FORS on the 
VLT we have obtained a spectrum of the brightest source 
N-15-2C, which lies at an angular separation of 1.53". This 
spectrum shows that the candidate source is a quasar at the 
same rcdshift as BR 2212— 1626 (the spectrum will be pre- 
sented in a future paper describing spectroscopy of our NIC- 
MOS and STIS candidates). While the hypothesis requires 
further investigation at present we consider it unlikely that 



this is a gravitational lens system, as we see no sign of a 
lensing galaxy in our NICMOS image. 

16. 2233.9+1318, Fig. 8. The faint diffuse source N- 
16-lD with Hab = 25.05 is the counterpart of the LLS at 
z = 3.1501. This object was first detected, and proposed 
as the couterpart, by Steidel, Pettini, and Hamilton (1995), 
who measured an offset from the quasar of 2.9", in agree- 
ment with otir value of 2.78". The object was discovered 
independently by Djorgovski et al (1996), who confirmed 
the identification through the spectroscopic detection of 
Lya emission. They quote a smaller offset of 2.3". This ob- 
ject has also been detected in our STIS image (M0ller et 
al, in preparation). As with the quasar H1500#13 after psf 
subtraction there is residual emission close to the quasar, 
inside our detection radius limit, which we do not believe is 
associated with the absorber. 



3.2.3 Comparison with other work 

Kulkarni et al. (2000) observed the quasar LBQS 1210+1731 
with the NIC2 camera and F160W filter for a single or- 
bit (2560 sec integration time), and detected a candi- 
date DLA absorber counterpart at an impact parameter 
of 0.25", at 4.1(7 significance. The object is compact and 
would have an aperture magnitude of Hab ~ 21.87. This 
object is rather bright. Placing it in Fig. 4 it can be seen 
that it is more than three magnitudes brighter than the 
three spectroscopcially confirmed DLA galaxy counterparts 
(Hab = 25.05, 25.11, 25.54), indicated by triangles. 

As seen in Fig. 4, our average detection limit becomes 
rapidly brighter at small impact parameters. Nevertheless 
their object lies well above our 6a detection limits. Our av- 
erage limit for compact objects at 0.3" is Hab = 23.87, with 
range 22.49 to 25.33 (Table 3), so wc would easily have de- 
tected a similar object in any of our frames. Indeed we do 
find a few bright candidates at small impact parameters. 
However we have chosen, conservatively, not to list candi- 
dates inside 0.56" for the reasons set out in section 2.4. 

The accuracy of psf subtraction of the two methods ap- 
pears to be comparable. We compute an equivalent detec- 
tion limit for their data of Hab = 22.12, by scaling their 
measured noise at 0.3" (0.28^Jy per pixel) to our longer in- 
tegration time. Their quasar has Hab = 16.95 ± 0.05 (total, 
P Hewett, private communication), similar in brightness to 
our two brightest targets, quasars 13 and 14. Our detection 
limits for these quasars are 22.72 and 22.49. 



3.3 Summary 

In this paper wc have described a search for the galaxy 
counterparts of 23 high-redshift high-column density Lya 
absorbers seen in the spectra of 16 quasars, using the HST 
NICMOS instrument. Within a box of side 7.5" centred on 
each quasar, over all the fields we have found a total of 41 
candidates, of which 3 have already been confirmed spectro- 
scopically as the counterparts. We provide detection limits 
as a function of impact parameter for each field, and the use 
of aperture magnitudes makes it very simple to compute 
whether a hypothetical galaxy of specified luminosity pro- 
file and impact parameter would have been detected in our 
survey. Our aperture-magnitude detection limits Hab ~ 25, 
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the small minimum impact parameter 0.56", and the sam- 
ple size make this the most sensitive search yet made for 
the galaxies producing high-redshift DLA absorption lines. 
At the same time we are obtaining very-deep optical images 
using STIS, which will provide an additional list of candi- 
dates (M0ller et al., in preparation). Except for unusually 
red galaxies the STIS images will reach substantially deeper 
than the NICMOS images, and therefore can be used to 
check the reliability of the candidate list provided in Table 
4. A discussion of the conclusions that may be drawn from 
the imaging data, and the limited spectroscopic follow-up so 
far completed, is deferred to the STIS paper. 

The next stage of this programme is a campaign of 
spectroscopy to measure the rcdshifts of the candidates, to 
identify which are counterparts of the absorbers. Since the 
galaxies are mostly faint, and close to the quasar, redshift 
determination will be difficult. The best hope will be to de- 
tect Lya emission. Although this lino is readily extinguished 
by dust, and therefore may be weak, it is likely to be the 
easiest to detect because at this wavelength the light from 
the quasar is removed by the absorber itself. 
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